Manipulating and protecting entanglement by means of spin environments 
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We study the dynamical behavior of two initially entangled qubits, each locally coupled to an environment 
embodied by an interacting spin chain. We consider energy-exchange qubit-environment couplings resulting 
in a rich and highly non trivial entanglement dynamics. We obtain exact results for the time-evolution of the 
concurrence between the two qubits and find that, by tuning the interaction parameters, one can freeze the 
dynamics of entanglement, therefore inhibiting its relaxation into the spin environments, as well as activate a 
sudden-death phenomenon. We also discuss the effects of an environmental quantum phase transition on the 
features of the two-qubit entanglement dynamics. 



I. INTRODUCTION 

The interplay between coherent and incoherent processes 
is key in the quantum mechanical processing of information. 
Systems designed in order to perform a given communication 
or computational task have to cope with the detrimental ef- 
fects of the surrounding world, which can affect an otherwise 
coherent process in many distinct ways 111]. In the context of 
distributed quantum information processing (QIP), where net- 
works of spatially remote quantum nodes are used in order to 
process information in a delocalized way, a good assumption 
is that each local processor is affected by its own environment. 
Such an architecture for a QIP device is currently at the focus 
of extensive and multifaceted theoretical and experimental en- 
deavors [Q]. 

Exstensive work has been performed on the incoherent dy- 
namics resulting from the coupling of QIP systems with baths 
weakly perturbed by the system-induced back-action. The 
loss of quantum correlations due to the environment has re- 
cently received considerable attention with particular em- 
phasis on the phenomenon of environment-induced entangle- 
ment sudden death |4l] (i.e., complete disentanglement in a 
finite time), which has also been experimentally tested for 
the case of electromagnetic environments |5]. And yet, es- 
pecially for solid-state implementations of quantum proces- 
sors, the case of structured environments is extremely rele- 
vant |@]. In this framework, the dynamics leading to complete 
disentanglement of two qubits coupled with a common spin 
environments (the so-called "central-qubit model") has been 
extensively studied ||7l-|lll- Exponential decay of the concur- 
rence fiol between two qubits initially prepared in pure states 
has been observed, the decay rate being enhanced for working 
points close to the quantum phase transition of the environ- 
ment LSfl. 

In this paper, we study the dynamical evolution of the en- 
tanglement between two remote qubits coupled with mutually 
independent spin environments. The exact time dependence 



of their reduced density matrix is obtained using an origi- 
nal approach lfl2i ITsIl which allows us to track the dynamics 
of quantum averages of one- and two-spin observables and 
entanglement properties. In stark contrast with most of the 
available literature, here we consider a "transversal" qubit- 
environment coupling that allows for energy transfer, result- 
ing in a dissipative-like behavior of the qubits and in a much 
richer entanglement dynamics. 

The transversal nature of the qubit-environment coupling 
allows us to reveal the occurrence of entanglement sudden 
death even when starting from initially pure states of the 
qubits, there included the maximally entangled ones which are 
extremely relevant for QIP. The above feature cannot emerge 
in the case of longitudinal couplings as reported in Refs. |i7|- 

We show that, by setting the environment in different op- 
erating regime, one can induce either entanglement sudden 
death or a freezing effect. Moreover, we shed light onto the 
differences in behavior experienced by the parallel and an- 
tiparallel concurrence introduced in |17], whose interplay is 
crucial in determining the properties of the quantum correla- 
tions within the two-qubit state, in particular at the environ- 
mental quantum critical point where the antiparallel entangle- 
ment appears to be better preserved than the parallel one. 

The paper is organized as follows. In Sec.|II]we define the 
physical system and the relevant interactions. In Sec. Hill we 
provide exact expressions for the dynamics of a single qubit 
coupled with a spin-environment. These results are used in 
Sec. |IV]in order to get the dynamics of the entanglement be- 
tween two qubits, each coupled with a spin-environment via 
a local isotropic (Subsection IIV Ab or anisotropic (Subsec- 
tion |IVBJ exchange interaction. Finally, in Sec.|V]we draw 
our conclusions and suggest a possible scenario where the 
main features of our physical model can be embodied. 
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FIG. 1 : Sketch of the physical model. Each of a pair of entangled 
qubits, and Qg, is locally coupled with a spin chain, and 
Fg, via the Hamiltonians "Hq^ and "Hog, respectively. The dynam- 
ics within each chain is ruled by the intra-chain Hamiltonians, "Hr^ 
and "Hrg- The wavy line indicates initial entanglement. 



II. THE SYSTEM 

We consider two non-interacting subsystems, A and B, each 
consisting of a qubit Qi^ (k = A, B) coupled with a chain F^^ of 
Nk interacting S -1/2 particles. Whenever useful we will 
hereafter use the index k - A, B so as to generically refer to 
either of the two subsystems A or B. As usual, the qubits 
are described in terms of S - 1/2 spin operators, which we 
indicate as so^- Operator s„^ (n^^ = 1 , A^^^) corresponds to the 
spin located at site n of the chain F^^. Notice that although the 
above notation suggests the spin describing Qi^ to sit at site 
of the respective chain F^^, this is just a useful convention, but 
has no implication about the physical nature of Qa and Qb- 

The intra-chain interaction is of XF-Heisenberg type with 
local magnetic fields possibly applied in the z-direction 

N^-l N, 



where /z„^ is the field applied at site «^ and J^^'J are the cou- 
pling strengths of the intra-chain interactions. Each F^ is 
open-ended, while neither the /;)^''s nor the magnetic fields 
h„^ need to be uniform along the chains. Qubit is coupled 
with the first spin of its environment, embodied by F^, via an 
exchange interaction of strengths J^^ and can be subjected to 
a local magnetic field /lo, directed along the z-direction. The 
corresponding Hamiltonian reads 



'Ho,^-2(J^Slsl + Jl^^)-2ho^sl^ 



(2) 



The Hamiltonian of the total system k is thus given by '/^ - 
'Hq^ +'Hr,. In Fig.[r|we provide a sketch of the model consid- 
ered. As stated previously, although A and B do not interact 
directly, they experience joint dynamics due to the possibility 
of sharing initial entanglement. Depending on the choice of 
the local interaction parameters and magnetic fields in "Hi^ the 
efforts required to tackle the model greatly changes. In Sec- 
tion |III] we describe the approach used in order to achieve an 
exact solution for the dynamics of local observables. 



III. EXACT SINGLE-QUBIT DYNAMICS 

We first concentrate on the dynamics of one subsystem only 
(thus dropping the index k throughout this Section). In partic- 
ular, we are interested in the evolution of a given initial state 
of the qubit Q, as determined by its coupling with the chain F 
and under the influence of the local magnetic field. We resort 
to the Heisenberg picture, which has been recently shown to 
provide a convenient framework for the analysis of quantum 
many-body systems of interacting particles [ 12]. The key step 
is the use of the Campbell-Baker-Hausdorff' (CBH) formula in 
the management of the time-evolution operator of the system. 
For an operator O associated with an observable of a physical 
system with Hamiltonian 'H, the CBH formula reads (we set 
h - I throughout the paper) 



(3) 



In virtue of the algebra satisfied by the Pauli matrices, we find 
that upon application of Eq. the time evolution of the com- 
ponents of So reads 



1 

1 



(4) 



^5(0 = -i2s'oit)f^(t), 

where cr^ (a = x,y,z) are the Pauli operators for the spin 
at site n and P„ - 111=1 The time-dependent coefficients 
n,J(f) and A^(f) are the components of the (A^-H l)-dimensional 
vectors n'(f) and A '(f) defined by 



p=() 



(2/7)! 



p=0 



(5) 



(6) 



where T stands for transposition, the vector v has components 
v'i - <5,o [W] and the tri-diagonal matrix r has elements 



■ 2hi6ij 



(7) 



The coefficients n''(f) and A'(f) are obtained from Eqs. Q 
and (|6]l by replacing t with . As tt^ is real and symmet- 
ric, there is an orthogonal matrix U that diagonalizes it, so 
that {tt^Y - UA}pU^, with A the diagonal matrix whose 
elements Aij - Ajdij are the (positive) square roots of the 
eigenvalues of tt^ . Similarly, there is an orthogonal matrix 
V that diagonaHzes t^t and such that (t^t)p = VA^'^V^ 
with the same diagonal matrix A as above. As a consequence, 
after straightforward matrix algebra, one can sum up the time- 
dependent series in Eqs. ^ and (|6|l to get 



IP(t) = Vn(t)V^v, A''(f) = VEifiV^v, 



(8) 
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where fl(t) and H(f) are diagonal matrices with elements 
Q,y(f) - cos(Ait)6ij and SyXO = sin(/l,f)<JiV- The above equa- 
tions hold regardless of the local magnetic fields or the cou- 
plings y„''* and y^ '. Adopting the language of Refs. fT2i[T3ll . 
the components of n ''''(f) and A-*'-''(f) embody the fluxes of 
information from the qubit Q to the spin chain F. 

By means of Eqs. (|4]i one can determine the time evolution 
of the single-qubit density matrix 

Po(0= ^2+X<^>)>'^o> (9) 

a 

where 1 is the identity operator and ( • ) indicates the expecta- 
tion value over the initial state of the system. Once the diag- 
onalizations required for determining the vectors n^'^''*(f) and 
A ''^'^(f) are performed, we need to evaluate the expectation 
values entering Eq. Such a task can be performed within 
two different scenarios. In the first, F has a small number of 
spins, so that one can design the precise structure of its initial 
state. In the second, F consists of a large number of spins, 
which puts the analysis in the thermodynamic limit, where we 
can benefit of specific symmetry properties of "Kr- Here, we 
concentrate on the latter situation, which has been the subject 
of several recent papers, due to the fact that it can be used for 
describing a proper qubit-environment system. 

We assume that Q and F are initially uncorrected and set 
the former in an arbitrary single-qubit state po and the latter in 
an eigenstate \^'^} of 'Hr- The initial state of the total system 
will thus be p = po (gi |*I'r) (*I'rl- Under such conditions, the 
procedure described in the above section is most conveniently 
implemented as Eqs. (|4]i greatly simplify due to the properties 
of "Hy- In particular, the conservation rule ["H, 07,] = 0, 
which states parity invariance of the Hamiltonian, together 
with the property {&"a^„) - for a fi and n + m, imply 

1 ^ 

<*S«)= 2 Z^'^O) [mOiPn^n) + A-mPnCT^,}] (10) 
1 ^ 

<*"o(0)= 2 [n;i(f)<A,'^„) - Al(t){Pn^:,}] (11) 

n=0 

«={) 

1 

-5 Z [n«(f)nm(0 + Kmiit)] {P„^yP,n&-:K) (i2) 

n<m 

1 ^ 

-- 2 [n;;(f)n;(f) + Aj,(f)A;(f)] {P„^,P„,d-l&-l) ■ 

Using Eqs. (l9ll-(fT2ll one can finally evaluate the single-qubit 
density matrix po(0- 

IV. EXACT TWO-QUBIT DYNAMICS 

We now consider the complete system AUB. Under the 
assumption of non-interacting subsystems, the time propaga- 



tor generated by "K^ + "Hb factorizes as 'llAif) ® 'tlB(f) with 
'U^if) - exp[-/'7Tf;ff]. Despite the absence of interaction, A and 
B might still display dynamical correlations depending on the 
initial state of the total system. In fact, if AUB is prepared in 
an entangled state, its dynamical properties will depend on the 
structure of such initial state and the entanglement evolution 
will follow from the interactions ruling A and B separately. 
The values of the parameters entering "K^ might thus be con- 
sidered as knobs for the entanglement dynamics. 
We prepare the total system at time f = in 

p(0) = p2^2«(0) ®/'H0) ®/''(0), (13) 

so that we can use the results of Sec.|lII]and write the two- 
qubit density matrix at time t in terms of the time-evolved 
single-qubit one. In fact, using the results of Ref. [ISJ, we 
have 

= Z ^il'^f^Pp^^^^ K^A,B) (14) 

with pfj'it) the elements of the single-qubit density matrix at 
time t and K a tensor of time-dependent coefficients. The 
two-qubit state can then be written as 

pT"(o - Z Z Ai:u('Kj:(^ppf"(^^ (15) 

PaJa PB.ra 

where lower-case indices take values 1 and 2, while capital 
ones are defined according to L = Ia + Ib - Ui^odl = 1, .., 4 
(with L - I,J,P,R and I - i,j,p,r). By setting pr,(0) = 
iTr,) {*I'r, |, i-e. by preparing both the chains in any of the 
eigenstate of their respective Hamiltonians, we can explicitely 
evaluate Eqs. (fT0l)-(fT2b and hence Eq. By comparing 
the latter with Eq. ( fT4l i we finally determine the coefficients 
K'"j{t), thus fully specifying the dynamics of the two-qubit 
state p2'' 28(f). 

Our approach is fully general and can be used in a variety 
of different situations. Here we concentrate on the case where 
the initial state of the two qubits is one of the Bell states 

|0±) = ^(|OO)±|11», 1^^) = ^(|01>±|10)), (16) 

which we dub parallel {\4>±)) and antiparallel (|(A±)) Bell 
states 1 17]. In virtue of the symmetries of such states, 
the concurrence of the two-qubit state can be written as 
C{f)-2 max{0, C||(f), C|x(f)), where we have introduced 

Cnit) = IP23|- VPI1P44 , (17) 
Ctt(0 = Ipwl- VP22P33 (18) 

with C|t(0) = 1 forp^"'^''(0) = \(p±){4)±\ and Ct|(0) = 1 for 
pQAQsf^Qj^ _ |(^^)((/f^|. Here the notation pu is used, for the 
sake of clarity, to indicate the matrix elements pfj^''{t) de- 
fined in Eq. ( fTSl ). In what follows Ca{t) [Cp(f)] is the concur- 
rence corresponding to the case where the qubits are initially 
prepared in an antiparallel [parallel] Bell state. 
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As for the environments, we take two identical chains of 
an (equal) even number of spins N with homogeneous intra- 
chain couplings and field, i.e. J'^,'^-J and h„^=h. Both chains 
are prepared in the ground state of 'T/p,, which is found via 
Jordan- Wigner and Fourier transformations lHlil. Straight- 
forward calculations yield the relevant mean values entering 
Eqs. Q as 



cos[(kf + I))?,,] sin[kf§„] 



i V sin I?,, / 

= /p p "X-x^ _ ^n,kF+l<Pm,kr - <Pn,kF'Pm,kF+l 
2(cosi?„ - COS17„,) 



(19) 



where, for convenience of notation, we omit the index k. We 
have introduced ipj^k- y/2/(N + 1) sm(j&t), ^k - kn/iN + 1), 
k e [1,N] and the Fermi wave number kf is determined by 
the magnetic field fl9ll . Moreover, due to the absence of sym- 
metry breaking in |*Pr,)i it is i^n^} — 0. We consider finite 
but large enough to avoid finite-size effects to influence our 
results. For the range of parameters considered here, N = 50 
is found to fulfill such condition and it is then chosen to set 
the length of both chains. Finally, we take the same coupling 
strength between each qubit and its chain, that is = 7y . 



A. Isotropic coupling between the qubit and the chain 



We now consider an isotropic coupling between each qubit 
and its respective environment, defined by 7q - Jq - Jo in 
Eq. (|2| In the theory of open quantum systems, such coupling 
typically corresponds to a dissipative interaction treated in the 
rotating wave approximation |20]. For isotropic coupling, the 
total Hamiltonians "/^ have rotational invariance along the z- 



axis. This implies r 



and thus W^(t) = Ui(t) = n„(f). 



(20) 



A^(t) - An(t) = A„(f), which allows us to write 
<5S(/)) = ^(no(t){&t;} + A„(f)<<J^,)), 

<^;,(f)> = ^(no(/)<<^;) - Ao(f)<<^i)), 

1 ^ 

n=0 

~ 4 X [n„(f)n„(0 + A„(f)A„(f)]^„m . 

From the above expressions we see that single-qubit states 
initially directed along the z-axis of the Bloch-sphere main- 
tain such alignment regardless of the Hamiltonian parameters. 
On the other hand, initial "equatorial" states with (i^(0)) - 
evolve in time and remain on the equatorial plane only for 
zero overall magnetic field. The rotational invariance around 
the z-axis of the total Hamiltonians "Ha and 'Hb has relevant 
consequences also on the evolution of the entanglement, as 
shown in Section lTVBI 

In the fully homogeneous case, i.e. for /iq = h and Jq - J, 
it is Aq = -2{h - cos and the /'' component of the cor- 
responding eigenvector has the same form as tpj^g with -i- 1 



replaced by H- 2. Hereafter, time will be measured in units 
of 7"'. In the thermodynamic limit where — > oo, the sum- 
mations can be replaced by integrals yielding 



<^o"(f)) = 
(Slit)} = 



2f 

^"1(20 r 



[cos(2/if)<o-o) - sin(2/if)<o-o)] ' 



2f 



[cos(2/i0<o^,> + sin(2/z0<o-o)] , 



(21) 



and, for the z component 



n=0 



njn even 



2r2 

u+ju (n + l)(m + 1) 



n,m odd 



4f2 



, (n + l)(m + 1) 
4f^ 



Jn+i(2t)J„,+i(2t)g„„, (22) 



J'«+l(2f)J'm+l(2f)^ 



where J'„{x) are the Bessel functions. Long-time expansions 
show that the mean value of the single-spin x(y) component 
decays as TJ. If ho - h - 0, Eq. (|22] ) reduces to (s'^it)) = 
y(f)((T^) with 7(f) = J^2f)l2f-, yielding a scaling at long 
times. Eqs. ^2\\ and (l22b give an excellent approximation also 
for finite A^ (< 50) within a time range where finite-size effects 
have not yet occurred. As the latter are caused by reflection of 
propagating excitations at the boundary of the chain, we can 
neglect finite-size effects for times up to ~ A^. In fact, as the 
maximum one-excitation group velocity is 2, it takes at least 
a time A^ for an excitation to leave Q and travel back to it. 

Let us now analyze the evolution of the entanglement be- 
tween the two qubits. We first notice that in the present case 
of isotropic coupling between and and given that the 
exchange interaction along the chain is set to be of XX type, 
both Hamiltonians have rotational invariance around the 
z axis. This enforces a disjoint dynamics of the off'-diagonal 
matrix elements (p23 and p^) entering Eq. ( fTTj i and ( fTSl ). As 
a consequence, if Ci-x(O) = it is C|j^(0 = at any time [the 
same holds for C||(0]- Moreover, we find that Ca(t)>Cp{t) for 
any symmetric t, irrespective of the Hamiltonian parameters. 
This means that antiparallel entanglement is more resilient to 
the eff'ects of the spin environment. Using the analytical solu- 
tions given above we finally obtain the time evolution of the 
concurrence. For h=0 and in the fully omogeneous case of 
Jq-J, we find 



C(0 = 2max^0,r2(f)+r(f)-^ 



(23) 



from which, by using the definition of y(f) in terms of Bessel 
functions, we infer that at Jesd - 0.9037 ESD occurs. In fact, 
our results show that sudden death occurs also for Jo + J, ex- 
hibiting coupling-dependent characteristics. It is remarkable 
that disentanglement at finite time occurs, here, for pure initial 
states of the two qubits, a feature due to the specific form of 
qubit-environment coupling considered here. 

In the weak coupling regime Jo <k J, Fig. |2] shows that the 
concurrence relaxation-time grows as Jo decreases. On the 
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C{t) 




FIG. 2: Concurrence between Qa and Qb for isotropic and weak 
coupling, Jq/J = 0.5,0.4,0.2 and 0.125 (going from left-most to 
right-most curve). No magnetic field is applied. 



(a) (b) 




FIG. 4: (a) Time ^esd after which the concurrence between and 
Qb first exactly vanishes (entanglement sudden-death), as a function 
of the coupling Jg/J (double logarithmic scale). No magnetic field 
is applied, (b) Entanglement revival time against logjo[io/^]- No 
magnetic field is applied. 



Other hand, for strong coupHng Jq » J, the non-Markovian 
character of the environments becomes evident [as shown 
in Fig. [3) and entanglement revivals occur due to the finite 
memory-time of the spin chain. These revivals can be in- 
tuitively understood as the result of the strong coupling be- 
tween the two-qubit system and the spins at the first site of 
each chain. A large Jq gives rise to almost perfectly coherent 
interactions within such qubit-spin pair, only weakly damped 
by leakage into the rest of the chains. A revival time (the time 
after which the revival exhibits a maximum) of about tt/Jq is 
inferred from Fig. [3] for Jo/ J » 1. In Fig.|4](a) such findings 
are summarized in terms of the dependence of log[o[fESD] on 
logjo[/o/-/]- The quasi-linear trend shown there reveals that 
the growing rate of Iesd for 7o/7 < 1 is slightly larger than 
the decreasing rate at Jq/J > 1. In Fig.|4](b) we show the be- 
haviorof the revival time logio[frev] againstlogio[/o/-/], which 
also exhibits a quasi-linear trend. 

The presence of finite magnetic fields significantly changes 
the dynamics of the entanglement. If ho > and h - Q, i.e. 
the field is only applied to the qubits, we expect an effective 
decoupling of each qubit from the dynamics of its environ- 
ment, such that both sq, precess with a Larmor frequency that 
depends mostly on ho, though it is subjected to small quanti- 
tative coiTections due to the interaction with F^^ at rate Jo- In 
fact, in Fig. |5] we see that a larger amount of entanglement 
is mantained for considerably longer times as ho grows. Due 
to the above decoupling, as well as the condition ho^ - ho^, 
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FIG. 3: Concurrence between (2a and Qb for isotropic and strong 
coupling, Jq/J = 3.5,4,4.5 and 5 (solid, dashed, dot-dashed and 
dotted curve, respectively). No magnetic field is applied. We also 
show the case of Jo/ J = 1 



correlations between Qa and Qb are preserved and so is their 
concurrence. 

In Fig. |5] (b) the same effect is illustrated by the time- 
averaged concurrence Ca,p — i^/St) J^^ Ca,p(t')dt' (the average 
is calculated over a time window 6t that excludes the oscilla- 
tory transients observed in Fig.|5](a) for Jt < 10 and Jt > 45). 
The average entanglement grows with ho- We further notice 
that parallel and antiparallel concurrences are almost identical 
(with Ca > Cp as expected), though their difference vanishes 
only as ho — > oo. We have also found that for /jq, hog the 
phase relation between individual precessions is lost and no 
entanglement preservation is consequently observed. 

We now switch off the magnetic field on both qubits (that 
is, we take ho - 0) and apply a finite field h > on the en- 
vironments. In this case a particularly interesting effect is ob- 
served as h becomes larger than the saturation value h - J 
and the dynamics of both chains slows down. As a conse- 
quence, after the transient, the dynamics of the correlations 
between the two qubits is considerably suppressed, due to 
the difficulty of the qubits to exchange excitations with sat- 
urated environments. A long-time entanglement memory ef- 
fect results from this, which is evident in Fig. |6](a). There, 
we also notice a reduction of the wiggling, which further wit- 
nesses the freezing of the entanglement dynamics. It should 
be remarked that such effect is profoundly different from the 
decoupling mechanism highlighted previously, where Ca-Cp 
was a monotonic function of the magnetic field. Here, in fact, 
a peak occurs in the difference between time-averaged con- 




FIG. 5: (a) Ca(t) versus /, for isotropic and omogenous coupling 
{Jo=J) on the chains with h,=Q and ;io/y=0.8, 0.9, 1, 1.1, 1.2, 1.5,2 
(from bottom to top curve), (b) Average parallel and anti-parallel 
concurrences and their difference plotted against Hq with the same 
parameters as in panel (a). The lines joining the data points are sim- 
ply a guide to the eye. 
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currence components when h-J, as shown in Fig. (|6]l (b), re- 
veaHng a drastic change in the entanglement behavior at the 
onset of an environmental QPT ll2lll . Clearly, at the envi- 
ronmental critical point, the antiparallel entanglement is fa- 
vored against the parallel one, which is at the origin of the 
peculiar behavior observed in Fig. |6](b) for the dashed line. 
The drastic change in the behavior of the average concur- 
rence observed at h/J - 1 is unique of the mechanism dis- 
cussed here and, as already stressed, clearly distinguished 
from the freezing effect due to mismatched frequencies at 
each qubit-environment subsystem. For h>J, the effect is 
fully established and the average concurrence increases, while 
Cfl and Cp get closer to each other. Moreover, by defining 
Z(t) = (l/2f2) Zn^i (n + \fjl^^{2t) and using the exact ana- 
lytical expressions (valid for h > J) 
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FIG. 7: Concurrence between Qa and Qg for anisotropic coupling, 
■^iSi ~ -^Ob ~ The magnetic fields are set to zero everywhere but 
on the chain Fa, where the field is let to change within the saturation 
region, = 1, 2, 10 (bottom to top). Solid (Dashed) lines are for for 
Cf^ (C-ji). All quantities are dimensionless. 



Ca(0= 2max<^0,y(f)- 



C/f) = 2 max j 0, j(t) + y'it) + Z\t)--^, 



1 [7Ht)+ZHt)] 
16 2 

1 



[yHt)-ZHt)f 



(24) 



we see that when the environments are saturated (i.e. all the 
spins of the chains are aligned along the z axis) the concur- 
rence dynamics does not depend on the magnetic field. 



B. Anisotropic coupling 

We finally consider the case of anisotropic coupling 7^ 7^ 
between each qubit and its respective environment (the chain). 
Differently from the case of isotropic coupling studied in Sub- 
section lTV Al and as a direct consequence of the fact that the to- 
tal magnetization of A and B is not conserved, the off-diagonal 
elements p23 and of the two-qubit reduced density matrix 
are not dynamically disjoint. This implies the possibility for 
the concurrence of the two-qubit state to switch from the par- 
allel to the antiparallel type and viceversa. 

For the sake of clarity, we consider extremely anisotropic 
conditions, setting = 0. In the case of no magnetic field 
on both qubits, /iq = 0, a very simple expression for the con- 
currence is found, due to the fact that (ij, (f)) is a constant of 
motion. In particular, if the two qubits are initially prepared in 




a combination of the two antiparallel Bell state, their concur- 
rence will evolve as C-rx(f) = -C-r-r(f) - (t)!]^ (t). On the 
other hand, if parallel Bell states are used to build up the ini- 
tial entangled state, it is C|i-(f) = -Cuit) - Yll^(t)TilJt). As a 

consequence, if 11^^(0 = 11^^(0 (holding whenTi^ = 'J^b), the 

one can drive the two- 



two-qubit concurrence cannot switch between C^i and C 
On the contrary, if U^i^(t) + 



qubit system from parallel-type to antiparallel and viceversa. 
In fact, the switching between parallel and antiparallel entan- 
glement is observed when so, and sog are flipped, under the 
effect of the coupling with the first spin of their respective 
chain, at different frequencies (for instance when 7^ + 7^), 
or when the dynamics of one subsystem is slowed down with 
respect to the other (for instance due to the fact that the field 
on one of the two chains is larger than the saturation value, as 
seen in Sec. lIV Al l. However, as discussed in Ref. |3], a two- 
channel entanglement evolution has an upper bound given by 
the product of the one-channel dynamics. Therefore, for such 
an "entanglement switching" to occur, the Hamiltonian pa- 
rameters of subsystems A and B should be set so as to re- 
tain high entanglement values. By fixing 7^, while setting 
/iA » 7^ in order to slow down the entanglement relaxation 
in the corresponding channel, the efficiency of the switching 
increases. This is clearly seen in Fig.|2l On the other hand, we 
can decrease 7;^ so that channel B is far more responsible for 
the entanglement evolution. In this case too a very efficient 
switching mechanism is achieved, suggesting that the satura- 
tion region of the chain is associated to an effective decoupling 
of the qubit from its corresponding environment. Finally, we 
notice that, being the coefficients defined by Eqs. ([8]l regular 
oscillating functions of time, the concurrence can only van- 
ish on a countable number of points on the temporal axis and 
cannot remain null for finite intervals of time. Therefore, en- 
tanglement sudden death is not observed. 



FIG. 6: (a) Cn(/) versus t, for isotropic and omogenous cou- 
pling (Jq=J), no field applied on the qubits (/io,=0) and h/J = 
0.8,0.9, 1, 1.1, 1.5,2, and 5 (from bottom to top curve), (b) Aver- 
age parallel and anti-parallel concurrences and their difference plot- 
ted against h/J with the same parameters as in panel (a). The lines 
joining the data points are simply a guide to the eye. 



V. CONCLUSIONS 

We have analyzed the dynamics of an entangled qubit-pair 
connected to two structured environments composed of open- 
ended and finite interacting spin chains. The intra-chain inter- 
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action has been modeled by an XX Heisenberg-like Hamilto- 
nian, while the coupling between each qubit and its respective 
environment has been realized via an XY exchange term with 
the first spin of the chain. Application of uniform magnetic 
fields has been also considered. We have exactly determined 
the time-dependent two-qubit density matrix, starting from in- 
formation gathered on the single-qubit dynamics. 

We have then provided analytical solutions both for the case 
of finite even and in the thermodynamic limit, thus getting 
access to a full-comprehensive and general analysis of entan- 
glement evolution. Particular emphasis has been given to the 
relaxation-like dynamics implied by the specific type of cou- 
pling considered, which gives rise, under suitable conditions, 
to a sudden death of the entanglement that we have analyzed. 
Interestingly, we have unveiled the occurrence of ESD also 
when starting from initially pure two-qubit states, a peculiar- 
ity of our model that does not emerge under "longitudinal" 
qubit-environment couplings. 

By manipulating the transverse magnetic field on the ini- 
tially entangled qubits, we have shown the possibility of de- 
coupling their dynamics from that of their respective environ- 
ments, thereby allowing for a dynamical entanglement protec- 
tion. On the other hand, when the magnetic fields applied to 
both chains are larger than the saturation value, the dynam- 
ics of the environments slows down and entanglement sudden 
death is not observed. Interesting features are observed when 
the environments undergo a quantum phase transition, which 
in this case happens when the field applied to the chain gets 
the saturation value, in particular as far as the the behaviour 
of parallel and antiparallel concurrence is concerned. 

Our work provides the analytical characterization of the 
transverse (i.e. energy exchanging) coupling between a simple 



and yet interesting out-of-equilibrium system (the two qubits) 
and a non-trivial spin environments (the two chains). As spin 
models are now understood as effective descriptions of many 
different physical systems, our results hold the premises to 
find fertile application to a variety of cases. As a particularly 
interesting situation, one can think about the engineering of an 
effective spin environment by using unidimensional arrays of 
small-capacitance superconducting Josephson junctions [li], 
which show a sharp phase transition from Josephson-type be- 
havior to Cooper-pair tunneling Coulomb blockade analogous 
to that of an XY model. This implementation thus constitutes 
a nearly ideal test for our predictions, since the effective envi- 
ronmental parameters can be modified through the use of gate 
voltages and external magnetic fluxes. It would be very inter- 
esting to study the applicability and relevance of a study such 
as the one performed here to the investigation of the properties 
of intrinsicaUy open systems in quantum chemistry and biol- 
ogy exposed to finite-sized environments. In this context, it 
is particularly significant that the mathematical model used in 
order to describe the radical pair mechanism B23ll bears some 
analogies with the central-qubit model. 
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